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ABSTRACT 

In  practice,  thyristor  rectifiers  are  used  as  high 
current-low  voltage  power  supplies.  Under  certain  con- 
ditions such  rectifiers  introduce  subharmonic  resonances 

In  this  paper,  a  limit  cycling  voltage  regulator 
system  is  simulated  in  the  digital  computer.   The  ref- 
erence voltage  and  load  are  varied  and  the  system  gain 
is  changed  for  better  understanding  of .their  effects  on 
the  subharmonics .   A  Fourier  analysis  of  the  output  wave 
shape  then  determined  the  existence  and  amplitude  of 
subharmonic  components. 
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I.   LIMIT  CYCLING  SYSTEMS,  THEORY  OF  OPERATION 

A.  INTRODUCTION  TO  LIMIT  CYCLING  SYSTEM 

In  practice,  it  might  be  required  to  have  a  control 
system  which  will  generate  a  limit  cycle.  In  this  thesis 
the  goal  is  to  investigate  the  characteristics  of  a  limit 
cycling  voltage  regulator.  This  kind  of  voltage  regulator 
system  is  potentially  useful  as  a  low  cost  source  of  well 
regulated  direct  current. 

The  system  introduced  in  this  paper  has  linear  and 
nonlinear  elements.   One  of  the  tools  for  linearizing  non- 
linear systems  is  the  describing  function  method.   In  deal- 
ing with  limit  cycling  systems  the  ordinary  describing 
function  method  predicts  the  limit  cycle,  but  is  not  ade- 
quate for  a  complete  analysis.   The  Dual  input  describing 
function  is  a  useful  method  for  investigation  of  the  limit 
cycling  system  when  subharmonic  resonances  are  anticipated. 

B.  DUAL  INPUT  DESCRIBING  FUNCTION  (DIDF) 

West  (4),  introduced  a  new  mathematical  device  which 
is  called  the  "Dual  input  Describing  function"  (DIDF) .   It 
is  an  extension  of  the  conventional  describing  function 
which  linearizes  the  nonlinearity  when  forced  by  two  inputs 

The  Dual  input  Describing  function  is  a  very  powerful 
tool  for  the  investigation  of  limit  cycling  systems.   It 
will  be  used  for  prediction  of  subharmonic  resonance. 
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For  better  understanding  of  a  limit  cycling  system  a 
general  block  diagram  is  introduced  in  Figure  1-la.   In 
this  limit  cycling  control  system  N(A,w)  represents  the  non- 
linear element,  G(juO  the  linear  element.   The  DIDF  is 
based  on  the  error  input  to  the  nonlinear  element.   For  a 
well  designed  system  the  error  X(t),  is  equal  to 

X(t)  =  r(t)-C(t) 

and  would  be  small  so  that  any  limit  cycling  period  To,  can 
be  modeled  as  shown  on  Figure  1-lb.   In  this  model  the  sin- 
usoid amplitude  A  is  associated  with. one  input  and  the  D.C 
bias  amplitude  is  associated  with  the  second  input.   Then 
the  input  to  the  nonlinear  element  is 

X(t)  =  B+AsinG 

Now,  the  calculation  of  the  DIDF  for  the  nonlinear  element 
N(a,to),  is  accomplished  by  computing  the  ratio  of  the  out- 
put amplitude  spectrum  of  the  nonlinearity  to  its  input 
which  gives  two  DIDF. 


N  (A,B)  =  ^  e"J6 


and 


°         A 


NS(A,B)  =  |1 


The  first  nonlinearity  is  considered  as  limit  cycle  DIDF 
and  will  be  used  to  predict  the  system  limit  cycle.  The 
second  DIDF  is  considered  as  a  linear  gain  for  range  B/A 
and  leads  to  linear  system  description  due  to  input.   The 
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FIGURE  1.1a.   Basic  limit  cycling  system. 
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FIGURE  1.1b.   Representation  of  input-output  relation  of  a 
nonlinear  element. 
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mathematical  definition  of  both  DIDF's  follow  directly  from 
Fourier  analysis  and  they  can  be  written  as 


Nq(A,B)    =    Re[No(A,B)]+JIm[No(A,B)] 


277 


TTA 


jp    J      y(B+Asin9)sin6d( 


and 


-^j-    J      y(B+Asin6)cos9 


27T 

Ns(AjB)  =  ~2iB~  /  y(B+Asin0)d9 


dO 


(1-D 


(1-2) 


Equation  (1-1)  can  be  reduced  since  integration  of  the  co- 
sine term  of  N  (A,B)  in  this  interval  is  zero.   Then  the 
simplified  form  of  N  (A,B)  can  be  written  as 


277 


No(A,B)    =      2ttB      /      y(B+Asin9)sin 


9d9 


(1-3) 


In  the  limit  cycling  system  B/A  will  be  small  .and  an  incre- 
mental DIDF  can  be  derived  as  B->o 


and 


DIDF  A  lim  [N  (A,B)] 


DIDF  = 


B+o 


1     lim 
277    B->o 


277 


-|-  J      y(B+Asin 


6)d6 


(1-4) 
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C.     DIDF  FOR  IDEAL  RELAY 

Let  us  consider  an  ideal  relay  with  drive  levels  ±D 
Figure  1.1c.   There  the  input  to  the  nonlinearity  is 

X  =  B+Asinu)0t 

=  B+Asine 

The  ideal  relay  has  a  single  valued  nonlinearity,  there- 
fore, computation  of  the  limit  cycle  DIDF  is  done  by  usinj 
the  result  obtained  in  equation  (1-3)  and  applying  it  to 
Figure  1.1c.   Then  the  DIDF  can  be  written 


2TT 

0 


2" 

N    (A,B)    =   — ±-      f      y(B+Asin0)sin0dl 

O  TTA         \s 


7T+6  !  27T-6i 

-4—  /  DsinSdO    +         J  (-D)sinGde 

ttA    *s  .  x 


A  - 


+  J      DsinBdf 

21T-0! 


COS01 


and 


ttA 


_i 

i    =    sin        (B/A) 


4D 


N    (A,B)    =   -S-        / 


ttA  v    1    -    (B/A)2  (1-5) 


16 


The  DIDF,  associated  with  the  signal  is 


VA>B)  ■  -iw 


ir+8, 

J         Dd9 


7T  +  0i 


J  (-D)d   +     J     Ddi 


2TT-I 


2D 


ttB 
2D 


sin    l    (B/A) 


(1-6) 


If  B-*-o ,    forumla    (1-6)    can  be   written 


N    (A,B)    =    lim    [N    (A,B) 
S  B+o         s 


VA'B>    =   "f- 


(1-7) 


And  the  DIDF  for  two  sinusoids  input  is  found  as 


X   =   Asin(u)   t  +  6)+Bsinoo    t 
o  <~> 


in  which  oo   is  the  limit  cycle  frequency  and  co   is  the  sig- 
nal frequency.   Then  the  DIDP  will  be 


■.<*.»>  =  8 


N  (A,B)  =  — 

S    '        TT 


3 


(  > 

it 

R 

A; 

— 

1+ 


s  (IJ    +    bT  (IJ 


(1-8) 


(1-9) 


If  this  result  is  compared  with  the  approximate  DIDF,  it  can 
be  said  that  the  approximate  DIDF  which  is  derived  as  equa- 
tion (1-5),  (1-6)  are  valid  for  5  percent  accuracy  but  if 
the  following  conditions  are  satisfied: 
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A    3 


-8<1 


Wo    3  (1-10) 

Then  it  is  possible  to  have  2.5  percent  accuracy  with  re- 
spect to  the  exact  DIDF. 

Therefore,  for  the  ideal  relay  case  the  limit  cycle 
DIDF  are  always  valid. to  better  than  5  percent  accuracy  re- 
ferred to  the  exact  DIDF,  under  the  restriction  of  (1-10). 

D.     LIMIT  CYCLING  SYSTEM  CHARACTERIZATION 

The  system  here  considered  is  a  kind  of  system  which 
always  has  a  limit  cycle.   Most  limit  cycling  systems  in- 
clude a  pre-filter  like  that  shown  in  Figure  1.2.   The 
fundamental  purpose  of  this  filter  is  to  reshape  the  command 
signal  amplitude  to  the  input  such  that  at  the  input  of  the 
nonlinearity  inequalities  (1-10)  are  satisfied.   This  brings 
a  design  problem  of  the  pre-filter  which  will  not  be  con- 
sidered in  this  thesis. 

1 .   Frequency  Domain  Description 

To  investigate  sinusoidal  frequency  analysis  of 
limit  cycling  systems,  the  system  model  is  illustrated  in 
Figure  1.2.  r(t)  is  a  sinusoid  and  the  input  to  the  non- 
linearity  N  is 

X  =  Asin(w  t+6)+Bsinw  t 
o  s 

The  goal  here  is  to  investigate  the  interaction  between  two 

sinusoids,  through  the  nonlinearity  and  differents  inputs 
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are  considered  in  this  analysis.   Figure  1.3  shows  the  range 
of  possibilities  of  the  signal  and  limit  cycle.   For  Fig- 
ure 1.3a  the  output  through  the  nonlinearity  is  nearly  sin- 
usoidal at  the  limit  cycle  frequency  but  is  phase  shifted 
with  respect  to  the  limit  cycle.   It  can  be  said  that  sig- 
nal through  relay  is  zero  which  means  N  =0.   Also  it  may 
be  said  that  for  signal  frequency  greater  than  the  limit 
cycle  frequency  but  satisfying  the  amplitude  ratio  con- 
dition also  gives  N  =0. 

In  Figure  1.3b  the  amplitude  ratio  is  satisfied 
but  the  frequency  ratio  is  0.5.   In  this  case  the  signal 
through  the  nonlinearity  would  not  lose  identity  and  the 
amplitude  ratio  only  will  depend  on  the  phase  shift  of  the 
input  component  to  the  nonlinearity.   Then  for  the  case 
shown  in  Figure  1.3b,  it  can  be  concluded  that  N  =N  ( 6 ) . 
For  a  frequency  ratio  equal  to  0.33,  N   affects  the  input 
signal  amplitude  by  factor  of  16  percent  and  a  phase  of 
4  degrees.   For  a  frequency  ratio  greater  than  0.33,  (but 
amplitude  ratio  satisfied)  N   can  be  considered  independent 
of  relative  time  position  since  the  change  in  amplitude  and 
phase  ratio  are  never  greater  than  1  percent.   Therefore  as 
a  result  of  these  conclusions  the  open  loop  transfer  function 
can  be  written  as 

CO 

Geq(Jw)=Ns.Heq.G(juO    for    |  u>  |  <   ^ 

0) 

Geq(joj)=Y.Ns.Heq.G(ja))    for   ^  <|a>|<u>0 
Geq(juO=0  for   ojo<|u)|  (1_1;L) 
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x=a  +  b 


FIGURE    1.3a.      AR   =    0.33,    FR   =    1.00 

x=a+b 


FIGURE    1.3b.      AR   =    0.33,    FR 

x=a+b 


0.5 


FIGURE    1.3c.      AR   =    0.33,    FR   =    0.16 
FIGURE    1.3.      Range    of  possibilities    of    input 
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where  y   represents  a  complex  number  depending  on  the  rela- 
tive phase  of  these  signals . 

2 .   Pole-Zero  System  Characterization 

Frequency  domain  analysis  shows  that  in  the  range 
between  -u  /3<oo  <co  /3  rad/sec  along  the  imaginary  axis,  if 
the  amplitude  ratios  are  also  satisfied  the  system  behaves 
linearly.   In  the  ranges  oo  /3<w  <w   and  -<o  <oo  <co  /3  rad/sec 

O   "™  S    O  O   S  ■"  o 

the  signal  gain  of  the  nonlinearity  is  a  function  of  time, 
and  for  the  remaining  portion  of  the  joo  axis,  the  signal 
gain  is  always  zero.   Now,  consider  exponentials  for  which 
the  amplitude  ratio  is  satisfied,  (amplitude  ratio  refers 
to  ratio  of  peak  exponential  to  peak  limit  cycle  amplitude) 
it  may  be  noted  that  the  exponential  excited  output  of  a 
relay  is  relatively  independent  of  the  time  origin  of  the 
input  exponential  signal  in  the  existence  of  the  limit  cy- 
cle, if  the  signal  lasts  for  two  or  more  limit  cycle  periods 
Figure  1.4  shows  this  situation.   Then,  the  minimum  value 
of  the  input  time  constant  is, 


4x  ,   «  2 

min 


2tt 


Then,  the  pole-zero  configuration  of  the  system  may  be 
described.   Consider  the  complex  plane  "p"  to  be  divided 
into  three  regions,  shown  in  Figure  1.5* 

The  nonlinear  element  N  of  the  system  will  be 
considered  as  a  gain  element  and  N  is  derived  for  the  relay, 
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FIGURE    1.4.      Response    of   ideal    relay    in   the   presence    of 
limit    cycle. 
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FIGURE  1.5.   Root  locus  plot  for  G(p) 
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therefore,  it  will  not  contain  any  phase  shifting.   If  the 
open  loop  system  gain  is  adjusted  to  one-half  of  the  value 
of  the  limit  cycle  gain,  it  is  possible  to  automatically  ob- 
tain a  limit  cycling  system.   The  closed  loop  roots  of  the 
system  can  be  approximately  determined,  such  that,  if  the 
gain  at  the  crossover  of  the  imaginary  axis  N=N  which  is 
the  limit  cycle  gain,  and  at  gain  value  of  N=N  /2=N   one 
of  the  closed  loop  roots  will  be  close  to  the  intercept 
point  of  the  to  /3  circle  on  the  real  axis  which  is  shown  as 

P  .   P  -co  /3  and  the  other  complex  pair  of  roots  would  be 
a    a  o 

very  close  to  the  circumference  of  the  to   circle.   This 

o 

result  determines  the  dynamic  response  of  the  limit  cycling 
system. 

One  of  the  important  points  about  limit  cycling 
systems  is  from  the  stability  point  of  view:   if  the  sensi- 
tivity of  the  system  is  increased  then  the  limit  cycle  am- 
plitude increases  but  the  gain  N   decreases  and  always  the 
S.N   product  remains  constant.   Closed  loop  root  locations 
do  not  change,  only  the  limit  cycle  amplitude  will  be 
changed.   Also,  it  is  possible  to  have  instability  in  the 
system  but  may  be  the  system  response  does  not  change. 

E.     CLOSED  LOOP  PROPERTY 

Consider  a  limit  cycling  system  as  shown  in  Figure  1.2 
Then  the  open  loop  transfer  function  of  the  system  is 


Geq(p)  =  Ns.H1(p).H2(p).G(p) 


=  Ns.Heq.G(p)  (1-12) 
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The  DIDF  for  the  ideal  relay  were  found  in  equations  (1-5) 
and  (1-6)  and  they  may  be  approximated  as, 

N   -  * 

o    ttA 

N   .  2£ 
s    irA 

and  the  ratio  of  these  two, 

No    2 

then  N  may  be  written  in  equation  (1-12) 
s 


G   (joo)  =  N  .H   .G(jco) 
eqw  '  s   eq    J 


IN  .H   .G(jo)) 
=  2   o   eq   VJ 


in  order  to  sustain  the  limit  cycle 


No.Heq.G(jM)  =  -  1  (1_13) 

must  be  satisfied.   Then, 

Geq(jw)  =  "  2  (1-14) 

This  is  a  kind  of  proof  of  the  explanation  given  earlier 
for  open  loop  gain  of  limit  cycling  systems. 
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II.   PERFORMANCE  CHARACTERISTICS  OF  THE  LIMIT 
CYCLING  VOLTAGE  REGULATOR 


A.     SYSTEM  DESCRIPTION 

The  block  diagram  of  the  voltage  regulator  system  used 
in  this  thesis  is  shown  in  Figure  2.1.  The  complete  circuit 
diagram  is  shown  in  Appendix  A. 

The  firing  circuit,  instead  of  using  a  separate  cir- 
cuit for  each  thyristor  uses  a  relay  to  turn  all  thyristors 
on  or  off.   If  the  relay  is  on  a  positive  voltage  is  applied 
to  all  thyristor  gates  and  this  causes  conduction  in  the 
rectifier.   If  the  relay  is  off  then  the  rectifiers  will 
not  conduct.  . 

The  rectifier  unit  is  a  three  phase  full-wave  recti- 
fier made  of  a  combination  of  six  thyristors.   The  thyristor 
is  a  solid  state  device  which  is  equivalent  to  a  thyratron. 
One  thyristor  device  is  a  silicon  controlled  rectifier.   It 
is  a  four  layer  device  and  may  be  likened  to  a  PNP  transis- 
tor joined  to  an  NPN  transistor  and  sharing- a  common  col- 
lector-base junction.   The  thyristor  and  the  characteristic 
curve  of  the  thyristor  are  shown  in  Figure  2.2. 

As  a  load  unit,  a  second  order  filter  with  very  low 
damping  factor  is  used.  It  is  an  LCR  filter.  The  values 
of  L,C  are  constant  and  variation  of  the  damping  ratio  is 
obtained  by  varying  the  load  resistance. 

The  system  has  a  nonlinearity  due  to  the  relay  and 
the  thyristor  characteristics  provide  a  second  nonlinearity. 
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FIGURE  2.1.   Basic  block  diagram  of  the  limit  cycling 
voltage  regulator. 
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FIGURE    2.2.      Thyristor   characteristic. 
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These  two  nonlinearities  are  treated  as  one  nonlinearity  In 
analysis  of  the  system,  providing  the  simplified  block  dia- 
gram of  the  system  as  shown  in  Figure  2.3. 

B.     STEADY  STATE  ACCURACY 

It  is  important  to  investigate  the  system  response  to 
different  inputs.   For  this  purpose  the  voltage  regulator 
system  is  analyzed  for  a  load  resistance  R=0.8ft,  and  dif- 
ferent reference  voltages  applied  to  the  system.   The  im- 
portant feature  to  be  studied  is  the  output  characteristic 
of  the  system.   For  different  reference  voltages  the  out- 
puts obtained  were  subjected  to  Fourier  analysis.   Table  2-1 
shows  the  result  of  this  analysis.   In  Table  2-1  results 
taken  from  digital  computer  calculation  of  the  Fourier  se- 
ries and  subroutine  FORIT  used  for  this  analysis.   In  the 

Table  2-1  A   is  the  DC  value  of  the  output. 
o 

These  results  show  that  the  system  output  is  not  equal 

to  the  required  value.   Then  the  difference  may  be  defined 

as  a, 

a  =  V   p  -  A 
ref    o 

For  high  reference  voltages  the  system  has  large  a  and  if 

V   f  gradually  decreases,  a  also  decreases.   At  V    =  52 

volt  a  changes  it's  sign,  that  is,  for  low  voltage  cases 

the  system  output  voltage  is  a  little  bit  higher  than  V   f • 

This  may  be  good  but  when  V   f  decreases  the  system  ripple 

instability  increases  which  is  not  desired  for  the  system. 

The  numerical  results  show  that  if  V   f  is  decreased  gradually, 
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FIGURE  2.3.   Simplified  block  diagram  of  the  limit  cycling 
voltage  regulator. 
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FOURIER      ANALYSIS     FOR      R=0.8fi 


ft£FERENCE 
VOLTAGE 
Vref(volT) 

A>it) 

UMiT  CYCLE 

flMPLiTUDE 

(volt) 

SUBHAMDNIC 

flMPLiTUDE 

(volt) 

UMiT  CYCLE 

TfttQUENCY 

(Hz) 

SlBHAiWONiC 

rftQUlNtt 

(Hz) 

ORDER    OF 
SllBHAWONIC 

64 

60.589 

3.34 

--- 

36.4 

--- 

1/2 

62 

59.685 

2.82 

0.373 

36.4 

18.2 

1/2 

60 

59.111 

2.36 

--- 

40.0 

--- 

1/2 

58 

57.24 

2.07 

0.146 

51.6 

25.8 

1/2 

56 

56.277 

2.4 

0.121 

51.2 

25.6 

1/2 

54 

53.961 

1.53 

0.134 

72.4 

36.2 

1/2 

52 

52.966 

2.21 

0.133 

60.4 

30.2 

1/2 

50 

50.272 

1.24 

0.143 

90.0 

45.0 

1/2 

48 

49.181 

1.505 

0.475 

56.2 

28.1 

1/2 

46 

47.497 

1.41 

0.635 

68.4 

17.1 

1/4 

TABLE       2-1 
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subharmonics  become  more  effective  and  this  may  make  the  sys- 
tem unstable.   For  V   f  =  58  volt,  the  amplitude  of  the  sub- 
harmonic  is  quite  small  (about  0.146  volt)  but  for  V   f  = 
46  volt,  the  subharmonic  amplitude  becomes  O.635  volt  which 
is  quite  large  compared  with  0.146  volt.   Also  the  results 
show  that  the  limit  cycle  amplitude  shows  an  opposite  re- 
action,  increase  in  the  subharmonics  amplitudes  caused  a 
decrease  in  limit  cycle  amplitude.   Figure  2.4  shows  the 
output  of  the  voltage  regulator  at  V   f  =  60  volt  and 

V   p  =  58  volt. 
ref 

C.     TRANSIENT  RESPONSE  OF  THE  SYSTEM 

For  investigation  of  the  transient  response  of  the  sys- 
tem only  the  linear  parts  are  taken  Into  account.   The  be- 
havior of  the  nonlinear  part  is  not  considered.   The  linear 
filter  of  the  system  is  assumed  to  be  as  shown  in  Figure  2.5 
And  the  transfer  function  of  the  system  is 

G(S)  -  Vout  1/LC 


Vin   "  S*+(1/RC)S+(1/LC)     (2-1) 

In  general  second  order  system  characteristic  equations  can 
be  written  in  the  form 


S2  +  2£oonS+ton2=  0 


then  for  the  filter, 


1 


and     ^  =  or  ^L 

/~tc~  2R   C 
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FIGURE  2.5.   Filter  and  load  of  the  system 
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During  analysis  of  the  system  L  and  C  are  held  constant  and 
only  the  value  of  R  is  changed 

C  =  O.lf. 
L  =  0.0005h. 
therefore, 

a)   =  141  rad/sec 

and  for  R  =  0 . 8ft 

i   =  0.044194 

which  is  a  very  low  damping  factor. 

The  settling  time  of  the  filter  is 

5     ^n 


Then, 


t   =  2.5    sec  for  R  =  3.2ft 
t   =  1.23   sec  for  R  =  1.6ft 
t   =  0.645  sec  for  R  =  0.8ft 
The  filter  transfer  function  can  be  written  in  Bode  form, 

G(jw)  =  LC(ja))2+L/R(ju3)  +  l 
And  for  R  =  0.8ft 

G(ja))  =   (1/2  10lt)(jco)2  +  (12.3/2  lO^CjuO+l 

(2-2) 
The  phase  and  amplitude  curves  of  the  Bode  diagram  are 
drawn  with  the  IBM-360  digital  computer  and  the  curves  are 
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illustrated  in  Figure  2.6.   The  analysis  of  the  closed  loop 
characteristic  of  the  filter  is  shown  on  the  Nichols  plot 
in  Figure  2.7- 

D.  PREDICTION  OF  LIMIT  CYCLES 

For  prediction  of  the  limit  cycle  frequency,  the  basic 

describing  function  method  is  used.   The  describing  function 

for  the  nonlinearity  of  the  system  which  includes  relay  and 

rectifier  is  derived  by  Leszcynski  (Ref.  15).   And  for  the 

linear  part  of  the  system  for  this  prediction  equation  -(2-2) 

is  used.   Figure  2.8  and  Figure  2.9  shows  the  Nichols  plot 

of  the  limit  cycling  voltage  regulator  for  V   f  =  64  volt 

and  Vref  =  62  volt  both  at  R  =  0.8ft.   Figure  2.8  gives  the 

limit  cycle  frequency  f,   =  35HZ  and  limit  cycle  amplitude 

B  =  3.6544  volt.   This  agrees  with  the  experimental  result 

shown  in  Table  2-1.   Also  for  V   „  =  62  volt  the  Nichols 

ref 

chart  of  Figure  2.9  gave  the  limit  cycle  frequency  f-.  =37HZ 
and  limit  cycle  amplitude  about  B  =  2.4632  volt.   Experi- 
mental results  in  Table  2-1  gives  f-,   =  36.4HZ  and  limit 
cycle  amplitude  about  B  =  2.82  volt  which  agree  fairly  well 
with  the  result  obtained  by  Nichols  plot. 

E.  SUBHARMONICS 

One  of  the  purposes  of  this  paper  is  also  to  investi- 
gate subharmonic  resonances  in  the  system.   Table  2-1  for 
Fourier  analysis  shows  relationships  between  V   ~  and  the 
subharmonic  amplitude  and  frequency.   Results  show  that 
changes  in  V   f  affect  subharmonic  amplitudes.   If  V   f  is 
increased  the  subharmonic  amplitudes  also  increase  and  for 
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FIGURE  2.6.   Bode 
diagram  of  the 
linear  part. 
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FIGURE  2.7.   Nichols  plot  of  the  filter  for  R  =  0.8ft 
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FIGURE    2.8.      Nichols   plot    of  system   for  Vref=    64V.,    R=0.8ft,    flc=35 
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GURE  2.9.   Nichols  plot  of  the  system  for  Vre£=62V.,  R=0 . 8Q 
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V   f  <  50  volt  subharmonics  become  very  important.   System 
behavior  gives  a  ripple  instability  at  the  output.   This  in- 
stability may  be  caused  by  the  high  amplitude  subharmonics. 
Figure  2.10  shows  this  kind  of  output  for  V   f  =  48  and 

Vr.of  =  k6     VOlt' 

rel 

Prom  the  frequency  response  point  of  view,  it  may  be 

concluded  that  changes  in  V   f  do  affect  appreciably  the 
frequency  of  the  subharmonics.   The  basic  frequency  of  the 
subharmonic  tends  to  remain  at  25  cps .   Table  2-1  shows  that 
changes  in  the  system  load  resistance;  R  =  0.8ft,  R  =  1.6ft, 
R  =  3.2ft  caused  no  change  in  the  subharmonic  frequency  and 
in  like  manner  changes  in  the  reference  voltage  did  not 
change  the  subharmonic  frequency.   These  results  shown  in 
Table  5-1,  Table  5-2,  and  Table  5-3.   It  may  be  concluded 
that  any  change  in  the  system  will  not  effect  the  subhar- 
monic frequency  since  the  system  forced  the  subharmonics 
toward  the  25  HZ  value.   But  the  amplitude  of  the  sub- 
harmonic is  quite  sensitive  to  any  change  in  the  system  or 
any  disturbance  to  the  system. 
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III.   INTRODUCTION  TO  THEORY  OF  SUBHARMONIC  RESONSNCE 

Subharmonic  resonance  is  a  phenomenon  which  occurs  in 
nonlinear  systems.   Usually  the  nonlinear  system  subharmonic 
resonance  frequency  is  a  fraction  of  the  input  frequency. 
As  an  example  consider  the  limit  cycling  voltage  regulator 
system  which  is  introduced  in  this  thesis,  or  a  relaxation 
oscillator,  or  a  moving-coil  loudspeaker  with  a  non-uniform 
magnetic  field  are  the  common  examples  of  the  subharmonic 
resonance.   This  phenomenon  can  also  occur  in  servo  mecha- 
nisms with  saturation  or  with  other  nonlinearities .   In  all 
of  these  cases  the  lowest  frequency  present  at  the  output 
is  a  submultiple  of  the  input  frequency  and  usually  at  or 
near  the  natural  frequency  of  the  system.   The  system  out- 
put amplitude  would  be  different  depending  on  the  existence 
of  the  subharmonic.   Such  subharmonics  can  increase  or  de- 
crease the  output  amplitude.   For  systems  in  which  the  out- 
put is  intended  to  reproduce  the  input  as  closely  as 
possible,  as  is  the  motion  of  the  loudspeaker  cone  or  a 
servo  motor,  this  phenomenon  can  introduce  an  intolerable 
distortion.   Also  it  is  observed  that  in  limit  cycling  vol- 
tage regulator  systems,  subharmonics  can  cause  instability 
of  the  system.   Also  one  of  the  important  points,  subhar- 
monic resonance  does  not  occur  in  pure  linear  systems.   If 
subharmonics  are  not  wanted  the  best  way  to  eliminate  them 
is  by  changing  the  parameters  of  the  linear  system  or  plac- 
ing a  compensation  element  rather  than  trying  to  change  the 
nonlinear  part  of  the  system. 
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It  is  possible  to  predict  subharmonic  resonance  in 
the  nonlinear  system  by  using  the  Dual  input  Describing 
function  and  by  solution  of  the  nonlinear  differential  equa- 
tion.  There  are  some  disadvantages  to  the  second  method. 
Since  the  output  of  the  nonlinear  element  is  expressed  as  a 
power  series  of  the  input,  successive  approximations  to  the 
correct  solution  converge  rapidly  only  if  the  amount  of 
nonlinearity  is  small.   And  also  the  solution  of  the  dif- 
ferential equation  does  not  give  a  physical  picture  of  the 
mechanism  of  subharmonic  resonance. 

Subharmonic  resonance  normally  exists  only  in  lightly 
damped  systems  with  nonlinear  energy  delays  or  with  nonlin- 
ear restoring  forces.   In  systems  giving  rise  to  subhar- 
monics,  the  nonlinear  element  is  frequently  characterized 
by  hysteresis.   Typical  nonlinearities  are  backlash,  mag- 
netic hysteresis  loops  and  relays,  etc.   It  is  also  impor- 
tant to  note  that  generation  of  subharmonics  is  dependent 
on  initial  conditions,  i.e.,  they  may  occur  in  a  certain 
range  of  amplitudes  and  frequencies.   If  in  a  system  a 
subharmonic  is  obtained  it  stays  "locked  in"  over  a  reason- 
able frequency  range.   If  the  driving  frequency  is  changed 
to  a  new  value,  the  system  output  frequency  changes  also 
to  a  new  value  which  is  co/n,  00  being  the  driving  frequency 
and  n  order  of  subharmonic. 
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IV.   DUAL  INPUT  DESCRIBING  FUNCTION  FOR 
PREDICTION  OF  SUBHARMONICS 


In  the  voltage  regulator  system,  power  supply  ripples 
feed  around  the  loop  and  interact  with  the  nonlinearity  of 
the  system  and  cause  instability  of  the  system.   The  limit 
cycle  resulting  from  this  instability  is  at  a  submultiple 
of  the  ripple  frequency  and  is  commonly  called  a  subhar- 
monic  resonance.   Normal  describing  function  theory  cannot 
predict  the  existence  of  the  subharmonic  resonance.   But 
by  using  the  Dual  input  describing  function  method,  it  is 
possible  to  find  a  mathematical  description  of  the  nonlin- 
ear part  of  the  system  which  can  be  used  to  predict  the 
existence  of  these  resonances.   In  general  the  forward 
transfer  function  of  the  voltage  regulator  system  can  be 
written 

Vout 


and 


=  NUw).G(ja)) 

(M-l) 


N(ja))  =  No(jw)+N1(joo) 


in  which  N  (jco)  is  the  DIDF  due  to  the  DC  input  to  the  non- 
linearity  of  the  system  and  N  (jw)  is  due  to  AC  input  to 
the  nonlinearity.   G(ju))  is  the  transfer  function  of  the 
linear  part  of  the  system.   If  E  represents  an  input  to 
the  nonlinearity  then  equation  (4-1)  can  be  rewritten, 

Vout  _  N  .G(j(i))+N1.G(ju) 

E  (4-2) 
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Frequency  analysis  is  accomplished  in  the  usual  fash- 
ion, by  applying  the  Nyquist  criterion.   The  result  could 
be  interpreted  as  an  encirclement  of  the  -1+jo  point. 
Therefore,  equation  (4-2)  can  be  drawn  in  vector  form  as 
shown  in  Figure  4.1.   It  is  seen  from  Figure  4.1  that  the 
N  .G(jco)  vector  is  placed  at  the  tip  of  the  N  .G(jw)  vec- 
tor.  Because  N  ( j co )  represents  the  DC  part  of  the  DIDF 
and  N1(jo))  represents  the  AC  part  of  the  DIDF,  N  (joj)  will 
have  one  component  but  N  ( j  co )  will  have  many  components  as 
the  phase  angle  changes.   But  during  this  variation  oj  =  to 

must  be  constant.   Then  if  (for  different  go  )  this  vector 

x 

plot  is  repeated  Figure  4.2  would  be  obtained.   The  N..  .G(ju>! 

vector  would  be  placed  on  the  tip  of  the  N  G(jto)  as  phase 

is  varied.   By  following  this  procedure  many  curves  can  be 

drawn  and  if  any  of  these  curves  N  .G(jco)  encircles  the 

-1+jo  point  then  in  Nyquist  sense  the  system  may  be  un- 
stable.  Such  a  system  may  have  a  subharmonic  resonance, 
but  this  is  not  guaranteed. 
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Na.  G(jcOx) 


FIGURE  4.1.   Vector  diagram  using  DIDF. 
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FIGURE    4.2.      Analysis    of   the   DIDF. 
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V.   ANALYSIS  OF  THE  EXPERIMENTAL  RESULTS 

A.     ANALYSIS  OF  LIMIT  CYCLING  VOLTAGE  REGULATOR 

The  goal  is  to  find  out  the  exact  cause  of  sub harmonics 
in  the  limit  cycling  voltage  regulator  system.   For  this  pur- 
pose the  IBM-36O  computer  was  used  for  solving  system  dif- 
ferential equations  and  subroutine  INTEG  2  is  used.   Appen- 
dix B  shows  this  computer  program  and  as  an  example  one  of 
the  outputs.   The  output  graph  was  obtained  for  three  dif- 
ferent loads  and  in  each  case  V   „  is  changed  from  64  volt 
to  46  volt.   Each  run  was  permitted  enough  time  so  that  the 
system  reached  steady  state.   From  these  output  plots  by 
analyzing  the  wave  shape  of  the  output,  the  results  in 
Table  5-1,  Table  5-2,  and  Table  5-3  were  obtained. 

Those  tables  showed  that  the  system  has  mostly  1/2 
order  subharmonics  and  a  few  1/4,  1/8  subharmonics .   Then 
it  may  be  said  that  the  limit  cycling  voltage  regulator 
system  has  an  even  order  nonlinearity .   Because  most  ex- 
perimental results  have  shown  that  nonlinearities  of  even 
order  give  even  order  subharmonics. 

In  the  Tables  amplitudes  are  given  peak  to  peak  in 
volts.   These  amplitudes  are  given  as  the  observed  values 
of  the  output's  highest  spike  and  lowest  spike.   These  low- 
est and  highest  spikes  become  smaller  if  the  load  resistance 
is  increased.   Increasing  load  resistance  caused  increase 
in  the  ripple  instability  at  the  output.   For  many  cases  it 
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ANALYSIS    OF  OUTPUTS    FOR     R=0.8  n 


REFERENCE 
VOLTAGE 

VftEF  (VOLT) 

LARGER 
AMPLITUDE 

UP- volt) 

SMALLER 
AMPLITUDE 

(PP.  volt) 

LIMIT  CYCLE 
FREQUENCY 

(Hi) 

SU3HAnMONK 
FREQUENCY 

(Hz) 

: ORDER  OF 
SUBHARMONU 

64 

8.66 

8.66 

29.8 

--- 

--- 

62 

7.1 

6.7 

35.8 

17.85 

1/2 

60 

5.9 

5.9 

39.7 

--- 

--- 

58 

4.72 

4.32 

50.8 

25.4 

1/2 

56 

5.5 

5.11 

50.8 

25.4 

1/2 

54 

3.5 

3.15 

68.6 

35.8 

1/2 

52 

4.72 

4.72 

57.8 

--- 

--- 

50 

2.76 

2.36 

87.8 

44.6 

1/2 

48 

5.31 

1.97 

56.2 

28.1 

1/2 

46 

11.8 

1.77 

68.4 

17.1 

1/4 

TABLE      5-1 
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ANALYSIS      OF     OUPUTS     FOR     R  =  l.6  n 


REFERENCE 
VOLTAGE 

VfttF  (von) 

LARGER 
AMPLITUDE 

(PP.  VOLT) 

SMALLER 
AMPLITUDE 

(PP.  volt) 

LIMIT  CYCLE 

FREQUENCY 

(Hz) 

SUBHARMONIC 
FREQUENCY 
(Hz) 

ORDER    OF 
SUBHARMONIC 

64 

5.3 

5.7 

27.8 

13.35 

1/2 

62 

3.54 

3.54 

35.8 

--- 

--- 

60 

1.97 

1.77 

57.8 

28.7 

1/2 

58 

1.18 

1.18 

134.0 

--- 

--- 

56 

0.985 

0.985 

115.5 

--- 

--- 

54 

2.76 

0.788 

115.0 

25.4 

1/4 

52 

3.15 

0.984 

145.6 

72.8 

1/2 

50 

3.15 

0.59 

57.8 

25.4 

1/2 

48 

3.55 

1.22 

73.68 

3.21 

1/8 

* 

46 

4.32 

1.37 

80.8 

40.4 

1/2 

*    Subharmonic    is    determined  by   Fourier    analysis 
TABLE      5-2 
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ANALYSIS     OF     OUTPUTS     FOR     R-3.2  n 


REFERENCE 
VOLTAGE 

LARGER 
AMPLITUDE 

(P?.VOlt) 

SMALLER 

AMPLITUDE 

(pp.volt) 

LIMIT  CYCLE 
FREQUENCY 

(Hi) 

SUBHARMONIC 
FREQUENCY 

(Hz) 

ORDER    OF 
SU3HARMONI 

64 

3.15 

3.15 

21.2 

--- 

--- 

62 

1.57 

1.48 

39.7 

19.85 

1/2 

60 

1.18 

0.788 

68.8 

36.2 

1/2 

58 

2.36 

0.394 

104.8 

52.4 

1/2 

56 

1.575 

0.985 

90.6 

25.5 

1/2 

* 

54 

1.77 

0.59 

100.0 

50.0 

1/2 

52 

1.575 

0.786 

91 

23 

1/4 

* 

50 

5.9 

0.394 

63.6 

17.4 

1/4 

* 

48 

2.75 

1.18 

62.88 

31.44 

1/2 

46 

3.94 

0.53 

49.28 

12.32 

1/4 

*Subharmonic   is    determined  by   Fourier    analysis. 
TABLE      5-3 
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was  not  possible  to  determine  the  order  of  the  subharmonics 
by  inspection  of  the  output  curves.   In  these  cases  a 
Fourier  analysis  was  made  to  establish  the  subharmonic  fre- 
quencies and  amplitudes . 

Figure  5- la  to  5.16.   shows  the  wave  form  at  the  filter 
element  and  from  Figure  5.2  to  Figure  5.34  shows  the  output 
and  input  plot  to  the  filter.   If  one  of  the  input  plots  is 
carefully  observed  one  notices  the  difference  in  width  of 
these  pulses.   The  system  has  sinusoidal  input  to  the  thy- 
ristors  and  the  output  of  this  rectifier  goes  to  the  fil- 
ter.  This  will  charge  the  capacitor  at  the  filter  and  when 
the  output  reaches  V   f  the  relay  would  be  open  but  at  this 
moment  the  rectifier  will  still  continue  to  conduct  and 
during  these  times  the  capacitor  will  charge.   Afterwards 
it  will  discharge  till  the  output  voltage  reaches  V   f. 
If  the  output  voltage  is  much  higher  than  V   f  then  the  ca- 
pacitor will  discharge  for  a  longer  time,  therefore,  the 
system  would  give  a  different  pulse  width  due  to  this  ef- 
fect.  This  can  be  seen  by  the  different  number  of  ripples 
at  the  top  of  the  input  voltage  wave  form. 

Output  plots  from  Figure  5.2  to  Figure  5.11  are  for 
load  resistance  R  =  3.2ft.   At  lower  voltages  deciding  the 
order  of  the  subharmonic  by  inspection  becomes  impossible. 
Those  outputs  were  analyzed  in  Fourier  point  of  view  and 
those  outputs  obtained  by  using  Fourier  analysis  are  desig- 
nated with  asterisks  in  the  Table. 
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FIGURE  5.1a.   Input  to  the  filter,  V  f   =  60,  R  =  1.6ft 
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0.1988    0.22      0.24      0.26      0.28     0.29813 

FIGURE  5.1b.   Inductance  current,  V  r  =  60V.,  R  =  1.6ft 

ret       * 
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FIGURE  5.1c.   Output  voltage,  Vref  =  60V.,  R  =  1.6ft 
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FIGURE  5.  Id.   Inductance  voltage  for  Vrg£  =  60V.,  R  =  1.652 
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Figure  5-12  to  Figure  5.1*1  shows  the  input  plot  to 
the  filter  with  average  value  of  the  pulse  width  approxi- 
mately corresponding  to  the  limit  cycle  frequency. 

Output  and  input  plot  for  the  R  =  1.6ft  case  are 
shown  in  Figure  5.15  "to  Figure  5-26.   For  R  =  0.8ft  case 
they  are  shown  in  Figure  5-27  to  Figure  5-3^. 

B.    ANALYSIS  OF  THE  D-C  SYSTEM 

In  this  paper,  one  of  the  important  things  to  be  deter- 
mined (if  possible)  is  the  exact  cause  of  the  subharmonics . 
Therefore,  to  get  results  it  is  necessary  to  make  some 
change  in  voltage  regulator  system.   Then  it  will  be  easy 
to  see  the  effect  of  the  element,  (i.e.,  relay,  thyristor) . 
For  this  purpose  the  rectifier  part  was  replaced  with  a 
D-C  source,  and  the  system  was  as  shown  in  Figure  5.35- 

For  analysis  of  the  system  the  C..5000  analog  com- 
puter was  used  and  the  system  output  obtained  for  R  =  0.8ft, 
R  =  1.6ft,  and  R  =  3.2ft.   In  this  system  the  filter  unit 
was  the  same  as  for  all  other  tests,  the  only  change  was 
replacing  the  rectifier  with  a  65  volt  battery.   Fig- 
ure 5.36  shows  the  analog  diagram  for  the  C.5000.   Output 
and  input  plots  for  different  loads  are  shown  in  Figure 
5.37  to  5.39.   These  outputs  show  that  the  system  did  not 
have  subharmonics.   Also  these  runs  were  made  for  different 
V   f  conditions  but  the  system  never  produced  subharmonics. 
This  means  that  if  the  system  has  a  nonlinear  element  but 
does  not  have  a  periodic  forced  input  then  the  system 
never  gives  subharmonics  at  the  output.   In  order  to  have 
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FIGURE  5.12.   Voltage  into  the  filter,  V   f  =  62V 
R  =  3.2ft 
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FIGURE  5.13.   Voltage  into  the  filter.   Vre£  =  56V.,  R  =  3.2ft 
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FIGURE  5.14.   Voltage  into  filter,  Vref  =  48V.,  R  =  3.2ft 
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FIGURE  5.25.   Voltage  into  filter,  Vre£  =  46V.,  R  =  1.6ft 
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FIGURE  5.26.   Voltage  into  filter,  Vref  =  60V.,  R  =  1.6fi 
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FIGURE  5.33.   Voltage  into  filter,  Vre£  =  60V.,  R  =  0 . 8ft 
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FIGURE    5.34.      Voltage    into    filter,    V      ,   =    50V.,    R   =    0.8ft 
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FIGURE  5.35.   Basic  block  diagram  of  D-C  case. 
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subharmonics,  a  control  system  must  satisfy  two  conditions, 
i.e.,  the  system  must  have  a  periodic  input  forcing  func- 
tion and  a  nonlinear  element.   From  this  it  might  be  said 
that  the  exact  cause  of  subharmonic  in  the  limit  cycling 
voltage  regulator  is  due  to  the  nonlinear  part  together 
with  relay  and  rectifier  and  also  the  periodic  input  to  the 
system  from  the  AC  ripple.   It  is  important  to  remember 
that  in  order  to  have  subharmonics  the  system  usually  must 
be  lightly  damped. 

The  block  diagram  shown  for  D-C  case  at  Figure  5.36's 
detail  of  the  elements  are  explained  in  Appendix  C.   The 
output  and  input  plot  of  system  shown  from  Figure  5-37  to 
Figure  5.39  are  for  steady  state  conditions.   Figure  5-^0 
shows  the  change  in  pulse  duration  of  the  relay  upon  the 
system  reaching  the  steady  state. 

C.    GAIN  CONSIDERATION 

As  a  last  investigation  of  the  limit  cycling  voltage 
regulator,  the  gain  was  changed  from  1  to  10,000  and  the 
output  plot  was  obtained  for  different  values  of  gain. 
Figure  5. 41,  Figure  5.42,  and  Figure  5-^3  show  these  out- 
put plots.   Output  shows  that  change  in  system  gain  did 
not  effect  the  output  wave  form  and  it  did  not  effect  the 
subharmonic  resonance. 
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FIGURE  5.41a.   Output  plot,  Vre£  =  60V.,  R  =  3.2ft,  Gain  =  5.0 
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FIGURE  4.41b.   Output  plot,  Vref  =  60V.,  R  =  1.6ft,  Gain  =  5.0 
7  of 


0.099461  0.12      0.14      0.16      0.18      0.2 
FIGURE  5.41c.   Output  plot,  Vre£  ■  60V.,  R  =  0.8ft,  Gain  =  5.0 
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FIGURE  5.42a.   Output  plot,  Vre£  =  60V.,  R  =  3.2Q,  Gain  =  1000.0 
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FIGURE  5.42b.   Output  plot,  Vref  =  60V,  R  =  1.6ft,  Gain  =  1000.0 
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FIGURE  5.42c.   Output  plot,  Vre£  =  60V,  R  =  0.8ft,  Gain  =  1000.0 
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FIGURE    5.43a.    Output   plot,    Vref   =    GOV,    R  =    3.2ft,    Gain   =    10,000 
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FIGURE    5.43b.    Output   plot,    Vre£   =    60V,    R   =    1.6ft,    Gain   =    10,000 
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FIGURE  5.43c.  Output  plot,  Vref  =  60V,  R  =  0.8ft,  Gain  =  10,000 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

Under  some  load  conditions  a  limit  cycling  voltage 
regulator  system  gives  subharmonic  resonances  at  the  out- 
put.  This  is  not  permissible  in  many  applications.   The 
goal  of  this  thesis  is  to  investigate  the  voltage  regula- 
tor system  and  to  try  to  find  out  the  causes  of  the  sub- 
harmonic  resonance. 

Removal  of  the  rectifier  and  substitution  of  a  bat- 
tery removed  two  basic  characteristics  from  the  system: 

1.  The  ripple  at  frequency  360Hz  due  to  the  recti- 
fication of  the  AC  supply. 

2.  The  delayed  shut-off  (i.e.,  continued  conduction 
after  the  error  detector  signal  gave  the  "off"  signal) 
due  to  the  thyristor  firing  characteristics. 

Since  the  subharmonics  are  observed  to  be  submulti- 
ples  of  the  limit  cycle  frequency,  therefore,  item  (1) 
cannot  be  the  cause  of  such  subharmonics.   Item  (2) 
occurs  at  the  limit  cycle  frequency,  and  results  showed 
that  the  input  voltage  to  the  filter  had  a  variable  pulse 
duration  for  rectified  AC  but  there  was  no  pulse  width 
modulation  which  probably  caused  the  subharmonics.   It 
can  be  seen  from  the  input  plots  that  the  ripple  begins 
at  different  voltage  levels  even  though  the  system  load 
resistance  and  reference  voltage  are  held  constant  during 


one  analysis.   The  thyristor  characteristics  apparently 
cause  this  pulse  width  modulation. 

Another  important  thing  is  the  prediction  of  sub- 
harmonic  resonance  by  using  the  Dual  input  Describing 
function.   If  the  Dual  input  Describing  function  is  derived 
then  by  designing  a  compensator  the  subharmonics  may  be 
removed  from  the  system.   As  a  further  study,  the  Dual 
input  Describing  function  of  the  nonlinearity  should  be 
derived,  taking  into  account  the  pulse  width  modulation. 
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INTER- PHASE 
TRANSFORMER 
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u Kj — O — \y 


THYRISTORS 


FILTER 


THREE  PHASE,  FULL  WAVE,  LIMIT  CYCLING  VOLTAGE  REGULATOR 
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APPENDIX  B:   DIGITAL  COMPUTER  PROGRAM  FOR  LIMIT  CYCLING 
VOLTAGE  REGULATOR 


The  constants  C(I)  represent 

C(l)  =  Maximum  value  of  input  sine  wave 

C(2)  =  Value  of  the  filter  inductance  is  henries 

C(3)  =  Value  of  the  filter  capacitance  in  farads 

C(4)  =  Value  of  load  resistance  in  ohms 

C(5)  =  Reference  voltage 

C(6)  =  Gain  of  amplifier  before  dead  zone 

C(7)  =  Magnitude  of  dead  zone 

C(8)  =  Gain  of  amplifier  after  dead  zone 
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EVEL   18  MAIN  DATE  =  70320  19/12/10 

DIMENSION  X(30) ,XD0T(30),C( 15) 
C(10)=1. 
IBP=1 

PI=3ol41592 
X(8)=IBP 
X(10)=OoO 
1  CALL  INTEG2(T,X,XD0T,C) 
X(ll  )=C(2)*XDOT(2) 
XMAX=C(1)*SIN( PI/3o ) 
THA=2«0*PI*60«0*T 
PHASEA=ABS(C(1 M  SIN(THA) ) 
PHASEB=ABS(C(1 ) v SI N< THA+2. 0*P I/3o ) ) 
PHASEC=ABS( C(l >+SIN(THA+PI/3o ) ) 
X(6)=AMAX1( PHASEA,  PHASE B,P HAS EC) 
IF(PHASEAoGEoXMAX) IPH=1 
IF(PHASEB.GEoXMAX) IPH=2 
IF(PHASEC.GEoXMAX) IPH=3 

206  ERR0R=X(1)-C(5) 
X(8)=IPH 
VK1=C(6)*ERR0R 
IF(VKlcGE.C(7) )VDZ=VK1-C(7) 
IF(VKloLToC(7) )VDZ=0o0 
IF(VKloLT.Cc0)VDZ=VKl 
VK2=VDZ*C(8 ) 

GO  T0(100,101,102,103,  104)  ,  IBP 

100  CONTINUE 

IFCVK2)    207,207,208 

207  IBP=1 
VIN=X(6) 
X(5)=VIN 
X(7)=0o 

X(3)=(1./C(2))*(VIN-X(D) 
XD0T(2)=X(3) 

IF(X(2)oLTo-0.01 )     X(2)=-0o01 
X«4)=(l./C(3))*X(2)-(lo/(C(3)*C(4) ) )*X(1) 
XD0T(1)=X(4) 

VLAST=VIN 
X(8)=IBP 
GO  TO  1 

208  IF(VINoLT.VLAST)  GO  TO  207 
GO  T0( 101,102, 103)  ,  IPH 

101  VIN=PHASEA 
IF(VK2)207,207,105 

105  CONT INUE 
X(7)=lo 
X(5J=VIN 

IF(VIN-o 75)300, 301, 301 

300  VIN=0o0 
X(5)=VIN 
IBP  =  5 

GO  TO  302 

301  IBP=2 

302  X(3)=(lc/C(2) )*(VIN-X( 1) ) 
IF(X(2).LTo-0.01)  X(3)=0o0 
XDOT(2)=X(3) 

.  IF(X(2)cLTo-0o01 )  X(2)=-0o01 

X(4)-<l„/C(3))*X(2)-<lo/(C(3)*C(4) ) )*X(1) 

XDOTd  )=X(4) 

VLAST=VIN 

X(8)=IBP 

GO  TO  1 

102  VIN=PHASEB 
IF(VK2)207,207,106 

106  CONTINUE 
X(7)=2o 
X(5)=VIN 

IF(VIN-o 75) 303,304,304 

303  VIN=0o0 
X(5)=VIN 
IBP=5 

GO    TO    305 

304  IBP=3 

305  X(3)=(l./C(2) )*(VIN-X(1) ) 
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EVEL   18 


MAIN 


DATE 


70320 


19/12/10 


IF(X(2).LTo-0.01 )  X(3)=0o0 

XD0T(2)=xm 

IF(X(2).LTo-0.0U  X(2)=-0o01 

X(4)=(lo/C(3))*X(2)-(lo/(C(3)*C(4)))*X(l) 

XDOTd  )=X<4) 

VLAST=VIN 

X(8)=IBP 

GO  TO  1 

103  VIN=PHASEC 
IF(VK2)207,207,107 

107    CONTINUE 
X(7)=3. 
X(5)=VIN 
IF(VIN-o 75)306, 307, 307 

306  VIN=OoO 
X(5)=VIN 
IBP  =  5 

GO    TO    308 

307  IBP  =  4 

308  X(3)=(1«/C(2))*(VIN-X(D) 
IF(X(2).LTo-0.01 )     X(3)=0o0 
XD0T(2)=X(3) 
IF(X(2)cLTo-0«01 )    X(2)=-0o01 

X (4) 'i 1. /C ( 3 ) ) *X ( 2 ) - ( lo / ( C (3 )*C< 4) ) )*X ( 1 ) 

XDOTd  )=X(  A) 

VLAST=VIN 

X(8)=IBP 

GO    TO    1 

104  VIN=OoO 
X(7)=4o 
X(5)=VIN 

X(3)=(l./C(2) )*<VIN-X< 1) ) 
IF(X(2).LTo-0„01 )     X(3)=0o0 
XD0T(2)=X(3) 

IF(X(2)oLTo-0. 01 )     X(2)=-0o01 
X(4)=(l./C<3))*X(2)-(lo/(C(3)*C<4)))*X(l) 
XD0T(1)=X(4) 

IBP=5 

IF(VK2 )309, 309,310 

309  IBP=1 

310  CONTINUE 
X(8)=IBP 
GO    TO    1 
END 
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GULER  M„Ao   FORCED  LIMIT  CYCLING  REGULATOR   N04 
ONE  RUN  IS  CALLED  FOR 

INPUT  DATA  RECORD 

ORDER  OF  EQUATIONS  =   2 
INITIAL  TIME         =  OoO 
FINAL  TIME  =  OolOOOE  00 

STEP  SIZE  =  0o2222E-04 

THE  NON-ZERO  CONSTANTS,  C(I),  ARE 

C(  1)  =  0o6500E  02 

C(  2)  =  Oo5000E-03 

C(  3)  =  OolOOOE  00 

C(  4)  =  OolbOOE  01 

C(  5)  =  0o6000E  02 

C(  6)  =  OolOOOE  01 

C(  8)  =  OolOOOE  01 

ALL  THE  INITIAL  CONDITIONS  ARE  ZERO 


THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 

TIME  X(  0) 

VIN  X(  5) 

VOUT         X(  1) 
CURRT        X(  2) 


NO  GRAPHS  ARE  REQUIRED 
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GULER    MoAo       FORCED    LIMIT    CYCLING    REGULATOR       N04 


TIME  VIN  VOUT  CURRT 

0:44440E-03  O^oilii    02  0o"ll435E    00  0°o52178E    02 

0:88880E-03  0o63848E    02  0.46759E    00  0o 10752E    03 

0ol3332E-02  0o64986E     02  0.10692E    01  Oo 16424E    03 

0ol7776E-02  0„64303E    02  0o19205E    01  0322052E    03 

0o22220E-02  0o61820E    02  0o30149E    01  0o27452E    03 

0o26664E-02  0o576C6E    02  0o43375E    01  Oo 32447E    03 

0o31108E-02  0.5O918E    02  0o58689E    01  0o37138E    03 

0o35552E-02  0o63281E    02  0„76088E    01  0o42028E    03 

Oo39996E-02  0ob4871E    02  0o95628E    01  0o46975E    03 

0o44439E-02  0o64645E    02  0o11729E    02  0o51799E    03 

0o48883E-02  0o626l0E    02  0ol4097E    02  0o56321t    03 

0o53326E-02  0o58821E    02  0o16649E    02  0o60365E    03 

0o57770E-02  0o588C5E    02  0»19361E    02  0o63884E    03 

0o62213E-02  0o626C0E    02  0o22221E    02  0o67445E    03 

0o66657E-02  0o64641E    02  0.25231E    02  0o71005E    03 

0o71101E-02  0„64873E    02  0o28388E    02  CQ74392E    03 

0o75544E-02  0o63289E    02  0o31680E    02  0o77432E    03 

0    79988E-02  0o59933E    02  0.35086E    02  0o79954E    03 

8:84431E-02  Oo57588E    02  Oo38581E    02  0o81826E    03 

0o88875E-02  O06I8C8E    02  0.42143E    02  0o83557E    03 

0n93^18F-02  0o642^7E    02  0o45773E    02  0o85267E    03 

0    97762E-02  0o64987E     02  0o4o^64E    02  Oo 86794E    03 

Ool0221E-01  0„63857E    02  0.53207E    02  0o87971E    03 

0cl0665E-01  0o60939E    02  0o56980E    02  »o88633E    03 

0    11109E-01  0o56315E    02  0o60757E    02  0o88623E    03 

011554E-01  Oo50115E    02  0o64507E    02  0o87797E    03 

0    119O8E-01  0o42511E    02  O068I88E    02  0o86025E    03 

0ol2442E-01  0.33717E    02  0.71758E    02  0.83200E    03 

0ol2887E-01  0^23980E    02  0o75167E    02  0o79239E    03 

0ol3331E-01  0ol3571E    02  Oo78365E    02  0o74087E    03 

Ool3775E-01  0o27820E    01  0o812^9E    02  0o67717E    03 

Ool4220E-01  OoO  0o83917E    02  0o60402E    03 

0    14664E-01  OoO  0o86198E     02  0o52839E    03 

Ool5108E-01  OoO  0^88132E    02  0o45089E    03 

0    15553E-01  OoO  0o89714E    02  0o37183E    03 

0ol5997E-01  OoO  0o90937E    02  0o29152E    03 

0ol6441E-01  OoO  0o91798E    02  0o21028E    03 

0       6886E-01  OoO  0.92295E    8?  frll&U    02 

0ol7330E-01  OoO  0o92427E    02  0o46328E    02 


Oo 

13571E  02 

Oo 

27820E  01 

Oc 

0 

Oc 

,0 

Oo 

,0 

Oc 

,0 

0, 

,0 

Oc 

,0 

0< 

,0 

0< 

>0 

0, 

,0 

0< 

,0 

0, 

,0 

0, 

,0 

0, 

,0 

0, 

>o 

0, 

,0 

0, 

,0 

0 

00 

0 

o0 

0ol7775F-01  OoO  0.92228E    02    -OolOOOOE-Ol 

0    18219E-01  OoO  0°?}?7?E    02    -0. 10000E-01 

0„18663E-01  OoO  0o91717E     02    -OolOOOOE-Ol 

0    19108E-01  OoO  0o91462E    02    -0.10000E-01 

0oll552E-01  Olo  0o91208E     02    -OolOOOOE-Ol 

0ol9996E-01  OoO  0o90955E    02    -OdOOOOE-01 

0°20441E-01  OoO  0o90703E    02    -OolOOOOE-Ol 

Oo20885E-01  OoO  0o90451E    02    -OolOOOOE-Ol 

Oo21329E-01  OoO  0o90200E    02    -OolOOOOE-Ol 

0o21774E-01  OoO  0»89949E    02    -OolOOOOE-Ol 
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GULER  MoA, 


FORCED  LIMIT  CYCLING  REGULATOR   N04 


TIME 


VIN 


VOUT 


CURRT 


0o22218E- 

01 

OoO 

0o22662E- 

■01 

OoO 

Oo23107E- 

■01 

OoO 

Oo23551E- 

-01 

OoO 

Oo23995E- 

-01 

OoO 

0o24440E- 

■01 

OoO 

0o24884E- 

-01 

OoO 

0o25329E- 

■01 

OoO 

Oo25773E- 

■01 

0.0 

0o26217E- 

-01 

OoO 

0o26662E- 

-01 

OoO 

0o27106E- 

-01 

OoO 

Oo27550E- 

-01 

OoO 

Oo27995E- 

-01 

0.0 

0c28439E- 

-01 

OoO 

Oo28883E- 

-01 

OoO 

0o29328E- 

-01 

OoO 

0o29772E- 

-01 

OoO 

0o30216E- 

-01 

OoO 

0o30661E- 

-01 

OoO 

Oo31105E- 

-01 

OcO 

0o31549E- 

-01 

OcO 

0o31994E- 

-01 

OoO 

0o32438E- 

-01 

OoO 

Oo32883E- 

-01 

OoO 

Oo  33327E- 

-01 

OoO 

0o33771E- 

-01 

OoO 

0o34216E- 

-01 

OoO 

0o34660E- 

-01 

OoO 

0o35104E- 

-01 

OoO 

0o35549E- 

-01 

OoO 

Oo35993E- 

-01 

OoO 

0o36437E- 

-01 

OoO 

0o36882E- 

-01 

OoO 

0o37326E- 

-01 

OoO 

0o37770E- 

-01 

OoO 

0o38215E- 

-01 

0.0 

0o38659E- 

-01 

OoO 

Oo3Q103E- 

-01 

OoO 

0o39548E- 

-01 

OoO 

0o39992E- 

-01 

OoO 

0o40436E- 

-01 

0.0 

0o40881E- 

-01 

OoO 

0o41325E- 

-01 

OoO 

0c41770E- 

-01 

OoO 

0o42214E- 

-01 

OoO 

0o4265PE- 

-01 

OoO 

0o43103E- 

-01 

OoO 

0o43547E- 

-01 

OoO 

0o43991E- 

-01 

OoO 

0o89700E 
0o89451E 


02 
02 


OolOOOOE-01 
Oo  10000E-01 


0«89202E  02    -OolOOOOE-01 

0.88955E  02    -Oe10000E-01 

0o88708E  02    -OolOOOOE-01 

0o88462E  02    -OolOOOOE-01 


0o88216E    02 
0o87971E    02 


Oo  10000E-01 
Oo  10000E-01 


0.87727E  02    -OolOOOOE-01 

0.87484E  02    -3„10000E-01 

0o87241E  02    -OolOOOOE-01 

0o86999E  02    -OolOOOOE-Ol 

0o86757E  02    -OolOOOOE-Ol 

0o86516E  02    -OolOOOOE-01 

0o86276E  02    -OolOOOOE-Ol 

0»86037E  02    -OolOOOOE-01 

0o85798E  02    -OolOOOOE-01 

0o85560E  02    -0o10000E-01 

0.85322E  02    -OolOOOOE-Ol 

0o85085E  02    -OolOOOOE-01 


0o84849E 

02 

-OolOOOOE- 

-01 

0o84613E 

02 

-OolOOOOE- 

-01 

0o84379E 

02 

-Oo  10000E- 

-01 

0o84144E 

02 

-OolOOOOE- 

-01 

0.83911E 

02 

-Oo  10000E- 

-01 

0.83678E 

02 

-Oo  10000E- 

-01 

0.83446E 

02 

-OolOOOOE- 

-01 

0o83214E 

02 

-Oo  10000E- 

-01 

0o82983E 

02 

-Oo  10000E- 

-01 

0*82753E 

02 

-OolOOOOE- 

-01 

0o82523E  02  -OolOOOOE-01 

0o82294E  02  -0o 10000E-01 

0o82065E  02  -0o 10000E-01 

0«81837E  02  -OolOOOOE-Ol 

O08I6IOE  02  -OolOOOOE-01 

0o81384E  02  -OolOOOOE-Ol 

0o81158E  02  -OolOOOOE-01 

0o80932E  02  -OolOOOOE-01 


0o80708E    02 
0o80484E    02 


Oo 10000E-01 
Oo  10000E-01 


0o80260E    02    -OolOOOOE-01 
0o80038E     02    -OolOOOOE-Ol 


0o79815E 

02 

-OolOOOOE- 

-01 

0o79594E 

02 

-OolOOOOE- 

-01 

0»79373E 

02 

-OolOOOOE- 

-01 

0o79152E 

02 

-Oo 10000E- 

-01 

0«78933E 

02 

-OolOOOOE- 

-01 

0o78714E 

02 

-Oo  10000E- 

-01 

0o78495E 

02 

-OolOOOOE- 

-01 

0o782  77E 

02 

-OolOOOOE- 

-01 
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GULER    McA, 


FORCED    LIMIT    CYCLING    REGULATOR       N04 


TIME 


VIN 


VOUT 


CURRT 


0o44436E- 

-01 

OoO 

0o44880E- 

-01 

OoO 

0o45324E- 

-01 

OoO 

0o45769E~- 

-01 

OoO 

0o46213E- 

-01 

OoO 

0o46657E- 

-01 

OoO 

0o47102E- 

-01 

OoO 

0o47546E- 

-01 

OoO 

0o47990E- 

-01 

OoO 

0o48435E- 

-01 

OoO 

0o48879E- 

-01 

OoO 

0o49324E- 

-01 

OoO 

0o49768E- 

-01 

OoO 

Oo50212E- 

-01 

OoO 

0o50657E- 

-01 

OoO 

Oo51101E- 

-01 

OoO 

Oo  51545E- 

-01 

OoO 

0o51990E- 

-01 

OeO 

0o52434E- 

-01 

OoO 

0o52878E- 

-01 

OoO 

0o53323E- 

-01 

OoO 

Oo53767E- 

-01 

OoO 

0o54211E- 

-01 

OoO 

Oo  54656E- 

-01 

OoO 

Oc55100E- 

-01 

OoO 

0o  55544E- 

-01 

OoO 

0o55989F- 

-01 

OoO 

0o56433E- 

-01 

OoO 

Oo56878E- 

-01 

OoO 

Oo57322E- 

-01 

0,0 

Oo  57766E- 

-01 

OoO 

Oo58211E- 

-01 

0*0 

Oo58655E- 

-01 

OoO 

0o59099E- 

-01 

OoO 

0o59544E- 

-01 

OoO 

0o59988E- 

-01 

OoO 

0o60432E- 

-01 

OoO 

Oo60877E- 

-01 

0.0 

0o61321E- 

-01 

OoO 

0o61765E- 

-01 

OoO 

Go62210E- 

-01 

OoO 

0o62654E- 

-01 

OoO 

0o63097E- 

-01 

OoO 

0o63541E- 

-01 

OoO 

0o63984E- 

-01 

OoO 

0o64428E- 

-01 

OoO 

0o64871E- 

-01 

OoO 

0o65315E- 

-01 

OoO 

0o65758E- 

-01 

OoO 

0o66202E- 

-01 

OoO 

0o78060E    02 
0.77843E    02 


Oo  10000E-01 
OolOOOOE-Ol 


0o77627E    02  -OolOOOOE-01 

0o77412F    02  -OolOOOOE-Ol 

0.77197E    02  -OolOOOOE-01 

0o76982E    02  -OolOOOOE-Ol 

0.76769E    02  -OolOOOOE-01 


Oo76555E  02 
0o76343E  02 
0o76131E  02 


0o75920E  02 
0o75709E  02 
0o75499E 


02 


0o10000E-01 
Oo  10000E-01 
OolOOOOE-Ol 

0o10000E-01 
OolOOOOE-Ol 
•OolOOOOE-01 


0o75289E  02    -OdOOOOE-Ol 

0o75080E  02    -OolOOOOE-01 

0o74872E  02    -OolOOOOE-Ol 

0o74664E  02    -OolOOOOE-01 


0o74456E 
0o74250E 


02 
02 


0„  10000E-01 
Oo  10000E-01 


0o74044E  02    -0„10000E-01 

0o73838E  02    -OolOOOOE-01 

0o73633E  02    -OolOOOOE-Ol 

0o73429E  02    -OolOOOOE-01 

0o73225E  02    -OolOOOOE-Ol 


0o73021E 

02 

-Oo  10000E- 

-01 

0o72819E 

02 

-OolOOOOE- 

-01 

0.72617E 

02 

-OolOOOOE- 

-01 

0„72415E 

02 

-0o  10000E- 

-01 

0o72214E 

02 

-OolOOOOE- 

-01 

0.72013E 

02 

-OolOOOOE- 

-01 

0o71813E 

02 

-OolOOOOE- 

-01 

0«71614E 

02 

-0„  10000E- 

-01 

0.71415E 

02 

-0o  10000E- 

-01 

0o71217E 

02 

-0o  10000E- 

-01 

0»71019E 

02 

-Oo  10000E- 

-01 

0o70822E 

02 

-Oo  10000E- 

-01 

0o70625E 

02 

-Oo  10000E- 

-01 

0e70429E 

02 

-OolOOOOE- 

-01 

0o70234E 

02 

-Oo  10000E- 

-01 

0o70039E  02    -0.10000E-01 

0o69844E  02    -OolOOOOE-01 

0.69650E  02    -OolOOOOE-Ol 

0.69457E  02    -OolOOOOE-01 


0o69264E  02 

0o69072E  02 

O06888OE  02 

0o68689E  02 

0o68498E  02 


Oo  10000E-01 
OolOOOOE-01 
Oo  10000E-01 
Oo  10000E-01 
OolOOOOE-01 


0o68308E    02    -OolOOOOE-Ol 
O068H8E    02    -OolOOOOE-01 
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GULER    MoA.       FORCED    LIMIT    CYCLING    REGULATOR       N04 


TIME  VIN  VOUT  CURRT 


0o66645E-01 

OoO 

0o67929E 

02 

-Oo 10000E-01 

0o67088E-01 

OoO 

0.67741E 

02 

-Oo 10000E-01 

0o67532E-01 

OoO 

0.67552E 

02 

-Oo 10000E-01 

0o67975E-01 

OoO 

0.67365E 

02 

-Oo 10000E-01 

0o68419E-01 

OoO 

0.67178E 

02 

-OolOOOOE-01 

0o68862E-01 

OoO 

0o66991E 

02 

-Oo 10000E-01 

Oo  69306E-01 

0o0 

0o66805E 

02 

-Oo 10000E-01 

0o69749E-01 

OoO 

0o66620E 

02 

-Oo 10000E-01 

0o70193E-01 

OoO 

0o66435E 

02 

-Oo 10000E-01 

0o70636E-01 

0.0 

0.66250E 

02 

-OolOOOOE-Ol 

0o71080E-01 

OoO 

O066O66E 

02 

-Oo 10000E-01 

Oo71523E-01 

OoO 

0o65883E 

02 

-Oo 10000E-01 

0o71966E-01 

OoO 

Oo65700E 

02 

-0o 10000E-01 

Oo  72410E-01 

OoO 

Oo65518E 

02 

-Oo 10000E-01 

Oo72853E-01 

OoO 

0o65336E 

02 

-OolOOOOE-Ol 

0o73297E-01 

OoO 

0.65154E 

02 

-OolOOOOE-01 

0o73740E-01 

OoO 

0«64973E 

02 

-OolOOOOE-Ol 

0o74184E-01 

0.0 

0.64793E 

02 

-Oo 10000E-01 

0o74627E-01 

OoO 

0o64613E 

02 

-Oo 10000E-01 

Oo75071E-01 

OoO 

0o64434E 

02 

-OolOOOOE-01 

0o75514E-01 

OoO 

0.64255E 

02 

-Oo 10000E-01 

Oo75958E-01 

OoO 

0.64076E 

02 

-Oo 10000E-01 

0o76401E-01 

OoO 

0o63898E 

02 

-Oo 10000E-01 

0o76845E-01 

OoO 

0.6372  IE 

02 

-Oo  10000E-01 

0o77288E-01 

OcO 

0.63544E 

02 

-Oo 10000E-01 

0o77731E-01 

OoO 

0.63368E 

02 

-Oo 10000E-01 

0o78175E-01 

OoO 

0.63192E 

02 

-Oo 10000E-01 

0o78618E-01 

OcO 

0.630I6E 

02 

-OolOOOOE-Ol 

0o790o2E-01 

OoO 

0o62841E 

02 

-Oo 10000E-01 

Oo79505E-01 

OoO 

0o62667E 

02 

-Oo 10000E-01 

0o79949E-01 

OoO 

0.62493E 

02 

-Oo 10000E-01 

Oo  80392E-01 

OoO 

0o62319E 

02 

-Oo 10000E-01 

0o80836E-01 

OoO 

0.62146E 

02 

-Oo 10000E-01 

Oo81279E-01 

OoO 

0.61974E 

02 

-Oo 10000E-01 

0o81723E-01 

OoO 

0.61802E 

02 

-OolOOOOE-01 

0o82166E-01 

OoO 

0.61630E 

02 

-Oo 10000E-01 

0o82609E-01 

OoO 

0.614595 

02 

-Oo 10000E-01 

0o83053E-01 

OoO 

0.61288E 

02 

-Oo 10000E-01 

0o83496E-01 

OoO 

O06III8E 

02 

-Oo 10000E-01 

0o83940E-01 

OoO 

0.60948E 

02 

-Oo 10000E-01 

0o84383E-01 

OoO 

0.60779E 

02 

-Oo lOCOOE-01 

0o84827E-01 

OoO 

0c60610E 

02 

-Oo 10000E-01 

0o85270E-01 

OoO 

0.60442E 

02 

-Oo  10000E-01 

0o85714E-01 

OoO 

0o60274E 

02 

-Oo 10000E-01 

0o86157E-01 

OoO 

0.60107E 

02 

-Oo 10000E-01 

0o8660lE-01 

Oo6l298E    02 

0.59940E 

02 

-Oo 10000E-01 

0o87044E-01 

0o64042E    02 

0.59776E 

02 

0ol9466E    01 

0o87488E-01  0.64999E  02  0.59628E  02  0o63648E  01 
0o87931E-01  0o64144E  02  0o59501E  02  0ol0951E  02 
0o88374E-01       0o61501E    02       0»59392E    02       0o14083E    02 
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GULER    MoA.       FORCED    LIMIT    CYCLING    REGULATOR       N04 


TIME  VIN  VOUT  CURRT 

O0888I8E-OI       0o57142E    02       0o59291E    02       0ol4l88E    02 
0o89261E-01       0o60285E    02       Oo59186E    02       0o13191E    02 


0o89705E- 

-01 

0.63489E 

02 

0.59085E 

02 

0ol5766E 

02 

0o90148E- 

-01 

0o64922E 

02 

0o59001E 

02 

0o20489E 

02 

0o90592E- 

-01 

0o64546E 

02 

0c58940E 

02 

0o25748E 

02 

0o91035E- 

-01 

0o62369E 

02 

Oe58901E 

02 

0o29914E 

02 

0o91479E- 

-01 

0o58454E 

02 

0o58874E 

02 

Oo  31394E 

02 

0o91922E- 

-01 

0o59155E 

02 

0o58847E 

02 

0o30281E 

02 

0o92366E- 

-01 

0o62813E 

02 

0.58822E 

02 

0o32318E 

02 

0o92809E- 

-01 

0o64719E 

02 

0o58811E 

02 

0o36851E 

02 

0o93252E- 

-01 

0o64821E 

02 

0e58823E 

02 

0o42278E 

02 

0o93696E- 

-01 

0c63116E 

02 

0o58859E 

02 

0o46969E 

02 

0o94139E- 

-01 

0c59650E 

02 

Oo58910E 

02 

0o49318E 

02 

0o94583E- 

-01 

0c57909E 

02 

0o58964E 

02 

0o48260E 

02 

0o95026E- 

-01 

0o62014E 

02 

0o59015E 

02 

0o49249E 

02 

0o95470E- 

-01 

0o64390E 

02 

0.5P078E 

02 

0o53074E 

02 

0o95913E- 

-01 

0c.64971E 

02 

0o59160E 

02 

Oo58153E 

02 

0o96357E- 

-01 

0c.63739E 

02 

0c59265E 

02 

Oo62858E 

02 

0o96800E- 

-01 

0o60731E 

02 

0o59387E 

02 

0o65574E 

02 

0o97244E- 

-01 

0o56550E 

02 

0.59513E 

02 

Oo  64756E 

02 

0o97687E- 

-01 

0o61094E 

02 

0.59632E 

02 

0o64211E 

02 

0o98131E- 

-01 

0„63936F 

02 

Oo  59756E 

02 

0o66849E 

02 

0o98574E- 

-01 

0o64994E 

02 

0.598  96E 

02 

0„71106E 

02 

0o99017E- 

-01 

0o64239E 

02 

Oo60055E 

02 

0o75366E 

02 

0o99461E- 

-01 

0o61694E 

02 

0o60230E 

02 

0o78007E 

02 

AT    4500     INTEGRATION    STEPS 
ONE    RUN    CALLED    FOR    HAS    BEEN    COMPLETED, 
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APPENDIX  C 

SUMMER  AMPLIFIER: 

Amplifiers  on  the  patch  board  are  addressed  as  an 
even  number  A000  to  A126,  and  are  represented  with  a  gain 
symbol  at  the  input  to  the  amplifier 


INVERTERS : 

Sign  changes  are  addressed  from  A200  to  A217  and  are 
represented  as 


INTEGRATOR : 

These  amplifiers  are  addressed  on  the  patch  board  as 
odd  numbers  (A001  to  A0^7)  and  are  integrators.  They  are 
represented  with  a  symbol 

IA007 


A007  and  A011  were  used. 
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COMPARATOR: 

A  kind  of  device  which  gives  Ov .  or  -6v.  logic  level 
according  to  the  comparison  of  the  inputs.   Represented 
with  symbol 


X 

Ca 

1 

n     \ 

Y 

Ca  N 
Ca 

/ 

such  that  if 


fCa  =  -6v- 

X+Y    >o 

— 

/  C  =  Ov. 
1    a 

fC  =  Ov. 
>     a 

X+Y<o     - 

— 

C          -6v. 

DIGITAL/ANALOG  SWITCH  (DA) : 

It  is  a  switch  and  is  represented  with  a  symbol 

J   LOGIC 
SIGNAL 


^OUTPUT 


The  switch  is  closed  if  logic  signal  is  -6v.  and  open  if 
logic  signal  is  Ov.   These  logic  levels  usually  come  from 
a  comparator.   Switching  time  is  approximately  equal  to 
1  usee. 
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